New Mechanism of Quantum Oscillations in the Superconducting Mixed State 

(H cl < B < # c2 ) 



00 



o 

IS 

■4— » 

CZ3 



C 

o 
o 



> 
in 

O 
O 

oo 

On 

-i— > 
o3 



o 
o 



X 



Lev P. Gor'kov 1 ' 2 and J.R. Schrieffer 1 
1 National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL 32310 
2 L.D. Landau Institute for Theoretical Physics, Russian Academy of Sciences, 117334 Moscow, Russia 

(Received ) 



We argue that inhomogeneity inherent to the presence of 
periodic supercurrents in the vortex lattice sorts excitations 
by energies into the ones that are spatially localized and those 
that perform motion along large Larmour orbits. This energy 
threshold results in a new mechanism for the de Haas-van 
Alphen oscillations which enhances oscillations at B <C H C 2, 
even for an isotropic superconductor with a constant gap. We 
suggest that the mechanism is of a general character and can 
cause the slow decay of the de Haas-van Alphen effect when 
the field, B, decreases below H C 2- 

PACS numbers: 74.20. Fg, 72.15. Gd, 74.60.-w 

The de Haas-van Alphen (dHvA) effect in the super- 
conducting (SC)-state was recently reported at magnetic 
fields, B, surprisingly smaller than the upper critical field 
value, H C 2 [1,2]. Theoretically, the dHvA-signal was ex- 
pected to decay faster (see discussion and references in 
[3]). In the experiments [1,2] the effect remains observ- 
able at B as low as 0.3 0.4iJ C 2, where the average SC- 
gap is thought to be too large for an electron to perform 
a circular motion along orbits with the large Larmour 
radius, rr, ~ vpj^c- 

Inhomogeneity of the current distribution and of the 
SC-order parameter itself cause serious difficulties for a 
rigorous theoretical treatment. Therefore in [3] an ex- 
treme limit of the lattice of isolated vortices, 



Co < d < 6 L ; H cl < B < H c2 



(1) 



has been chosen to find out whether any symmetry zeroes 
in the gap function could enhance the dHvA effect in the 
SC-state (in (1) £o is the coherence length, d- the vortex 
lattice period, 5l- the penetration depth). 

The results [3] confirmed expectations [4] that in the 
3D-case existence of a symmetry line of zeroes in the gap 
may restore the dHvA effect for some field directions, 
although the effect is much weaker than predicted in [4] 
due to scattering of electrons on flux lines. 

As for a "d-wave" gap in layered superconductors, the 
levels' systematics is such [3] that it forbids their cross- 
ing the chemical potential with the field change - the 
mechanism known to be the essence of the dHvA effect 
in normal metals. No qualitative difference is expected, 
hence, between a "d-wave" and an ordinary anisotropic 
gap as far as the dHvA effect is concerned. 

Chances to observe the dHvA oscillations for a super- 
conductor in the regime (1) are rather vague, for the sig- 
nal rapidly decreases below the experimental resolution. 



Common predictions for the signal's amplitude would 
give an exponential factor, exp(— A/lo c ) <C 1, where A is 
a gap scale and uj c the cyclotron frequency (e.g., see [5]). 

On the theory part however, the limit (1) suggests sig- 
nificant simplifications, because the spatial distribution 
of the field, currents, and the gap are all well known from 
the phenomenological consideration [6]. Main contribu- 
tions are expected to come from the "bulk" since the 
volume occupied by vortex cores is small. 

In that which follows the field range (1) is chosen again 
to demonstrate the existence of some new quantum oscil- 
lations mechanism which is specific to the SC-state and 
is originated by spatial variation of the local density of 
states. Although we expect that the effect should have 
a more general character, we have chosen the regime (1) 
and the isotropic gap for its simplicity. 

In the SC-state an excitation bears electron- or 
hole-like features to the extent its energy e(p) = 
\/ v 2 F (j> — pf) 2 + A 2 , exceeds A. Taking supercurrents 
into account, the local spectrum becomes 



E(p) =e(p)+p-v fl (r) 



(2) 



where v s (r) stands for distribution of superfluid velocity. 
Hence, the second term in (2) generates a potential relief, 
such that an excitation with a momentum p may, or may 
not, perform infinite motion depending on whether its en- 
ergy exceeds, or not, the value A + max|p- v s (r)|. An ex- 
citation (2) experiences a Lorentz force on the part of the 
magnetic field. Whether an excitation may perform an 
itinerant motion is equivalent to the condition that such 
excitations can encircle a large Larmour orbit. The addi- 
tional potential in (2), thus, imposes an energy boundary 
between "localized" and "extended" states. The latter 
may contribute to magnetization. 

To show that the above simple arguments, indeed, lead 
to a new contribution to quantum oscillations in the SC- 
state, we explore further the equations derived in [3] for a 
superconductor in a magnetic field. The method [3] con- 
sists of averaging of the Gor'kov equations over classical 
trajectories in the magnetic field. 

The core of the method is given by eqs. (24-27) of 
[3]. We use the same notations [3] below and start with 
deriving the new Schrodinger equation with the h(tp)- 
terms (eq. (31) [3])included and A =const.: 



-■■ -2Eh(ip)y 



(E 2 - A 2 )y 



(3) 



Note that the h((p)-term of eq. (31) [3] is nothing but 
the Doppler shift in (2). In the derivation of (3) terms 



1 



of the order of h 2 (ip) and ui c (dh{ip) / dip) were omitted. 
To estimate the second term recall that the dependence 
h(ip) on ip originates from the dependence of v s (r(t)) on 
t: t — > <p/u! c . An electron moving along the orbit with 
a velocity Vf senses changes in v s (r) with a frequency 
v F /d. A change on distances ~ d is therefore equivalent 
to its change on the scale, Sip, of the order of 

Sip ~ (duj c /v F ) ~ d/r L (4) 

The /i(<p)-term being h(ip) <~ wf/^, one obtains: 

fcV) ~ w^'fa) ~ i^/d 2 ~ A 2 (£ /d) 2 

In (3) at E <~ A one has: 

Eh(ip) <~ A 2 (£o/d) < A 2 (5) 

As it was argued in [3] , the oscillatory part of the mag- 
netization is contained (eq. (6') [3]) in the expression: 

M = -£- c J2[KW\n(E,) (6) 

where the index A numerates the eigenvalues. At 
T = n(E\) = 1 (all states below zero are occu- 
pied). The bar (...) means the normalization integral, 
{2w)~ 1 f Q n (. • -)dip. Oscillatory effects, if any, come from 
the vicinity of the chemical potential [3] : 

/i = u c N + p, (7) 

The chemical potential being large, fi/u3 c 3> 1, the oscil- 
latory part of magnetization in (6) does not depend on a 
specific N . The periodic (in B~ x ) pattern in M of eq. 
(6) may, therefore, be expressed as a function of 

K = p,/w c (8) 

with k varying in the interval (0,1). The eigen functions 
y(ip) themselves are not periodic. The periodicity condi- 
tions are to be imposed on: 

y(ip)exp{~iKp} (9) 

With (8,9) in mind, the problem of solving eq.(3) be- 
comes equivalent to the one of finding the band structure 
of a particle moving in the periodic potential, 2Eh(ip), 
where k (8) plays the role of a quasimomentum, lu 2 =^> 
(2m e //) -1 and E 2 - A 2 being the "energy". 

Summation over A in (6) is routinely replaced by the 
integration by making use of the Poisson formula: 

+ OC +OC 

S(X-n)= exp(27riif A) (10) 

n= — oc K= — oc 

Applying (10) to eq. (6), integrating over A by parts, one 
would arrive at terms of the form: 



first noticed in [3]. Eqs. (10) and (11) have no value 
unless a relation between A and energy is established. 
Below we will construct a function 

\(E) = <5>{E)/2tt (12) 

such that all the eigenvalues in (6) are determined as in 
(10), by the provision: 

X(E X ) = n (13) 

If $(E) were known on the real E-&xis, it then may be 
analytically continued onto the complex plane. 
With this in mind first simplify notation in (3) 

- e = E + A - h max (ip ) (14) 

hmax(ipo) is the maximum of h{ip) along a given trajec- 
tory. One may assume ipo = 0. Eq. (3) becomes: 

- u? c y" + 2A(h(ip) - h max )y = 2(-e)Ay (3') 

Even though h <~ vp/d <~ A(^ /d) <g; A, the ratio 

Ah/u 2 c > 1 (15) 

is large. In the WKB-approach: 

y(<p) = a v+ if) + b y- O) (16) 

where y±{ip) are of the form [7]: 



y± (ip) = A(S')~ 1/4 




Here A is the normalization factor and 

S'(ip) = (v^A/wc) y/hmax - Hp) - e (17) 

At (— e) > the quantization condition (taking (9) into 
account) gives: 

V2A /" 27r 

5" (27T, -s n ) = / V h max - h(ip) - s n 

Jo 

= 2irn + 2hk (18) 

Eq. (18) at |e| » A matches the spectrum of free elec- 
trons in the magnetic field (in presence of the flux cur- 
rents). At (— e) < one would obtain "localized" states. 
Let us justify first this last statement. 

The "potential," h(ip) — h max < in (17) is rather 
irregular and has many minima alternating with local 
maxima on distances of order of 5ip, (4). The probability 
to tunnel between minima is measured by the value of 

1/2 

exp(— \5S\) with 5S ~ (<v£o) S> 1, i- e - is small enough 
to consider barriers as impenetrable ones. 



The two WKB-branches, one at (— e) > and the other 
for (— e) < 0, cannot together form a function Q(E) such 
as eq. (12). The WKB-approach is invalid at (—e) too 
close to the maximum of h(ip). Indeed, expanding (18) 
at small (— e) > results in the singularity: 



S(2n,-s) ~ 5(271-, 0) + 



(-£) 
2w c a 1 /2 



In 



Aa(6<pf 



(19) 



For h(ip) close to the maximum we have chosen: 

h(<p) = h max - (aAp> 2 /2) (20) 

with aA ~ vp / d{5tp) 2 . 

Thus, in the WKB-approximation there are two sorts 
of solutions: those of the form of eq. (16) at (— e) > 0, 
and the "localized" states at (— e) < 0. It is useful to 
show how these findings arc connected with our intuitive 
expectations regarding the role of the Doppler term (2) . 
The solutions u\(ip) and v\(tp) must be normalized 

1. There are two useful 



together: |ua(v)| 2 - 
auxiliary relations: 



H = (1/2) {|y|2 + (iuj c /2E)(y*y' - yy*')} 

W\ = (1/2) {W - {iuJ c /2E){y*yi -yy*')) (21) 

which may be derived with the same accuracy as (3). One 
concludes from (21) that \y\ 2 = 1, and |ua| 2 = \ inde- 
pendently on the energy for all localized states for which 
the wave functions are real. Therefore there is, indeed, a 
threshold singularity in M osc (11)- "localized" states give 
no contribution into M osc . To study the phenomenon rig- 
orously one must go beyond the WKB-accuracy. 

At the energy close to the maximum of the potential 
in cq. (3') one may apply the WKB solutions (16, 16') 
only far away from points ip = and ip — 27r. In the 
vicinity of (p — cq. (3'), (20) is solvable in terms of the 
parabolic cylinder functions. Matching the asymptotic 
behavior of the latter to the WKB-form (16, 16') makes 
it possible to find relations between coefficients (a, b) in 
(16) to the right of p = 0, and the similar set (a', b'), to 
the left of it: 



a [3 
(3* a* 



) (;) ; (N 2 -|/3| 2 = i) (22) 



Moving along ip from <p = (0)+ and using (16), one 
reaches close to the point ip — (2ir)-. With the help 
of (22), the solution then may be transformed into the 
solution at (p = (2ir) + . 

The periodicity (see (9)) leads to the equation: 

R(l) = \a\ (e +iS + e _iS ) =2cos27r« (23) 

Here and below the dimensionless energy, /, is 

I = 2(-e)/uj c a 1 / 2 (24) 



In (23) we used notations: 

a = \a\ exp(i6) ; S = S(2n,l) - 6(1) 



(25) 



(S(2ir, I) coincides with (19) re-written in terms oil from 
(24)). Of the two solutions in eq. (23) we have to take 



cos 2-kk, 



+ 



a 



'C0S27TK 



a 



(26) 



Indeed, at I — > +oo, a => 1, eq. (25) for energy levels 
goes over into the WKB-result (18). 
Introduce the function 

*(-e) = $(0 = S(2n, I) - 6(1) - i \np(l) (27) 

As we just mentioned, at large I solutions of $(/„) = 2irn 
go over into the WKB-spectrum which at very large ener- 
gies, / 3> 1, transforms into the spectrum of free electrons 
in the magnetic field. $(/) is a "monotonous" real func- 
tion at I ~ 1. At I < eq. (13) also determines the 
WKB-localized levels immediately below h max . There- 
fore we may use (27) as the definition of X(E) (12), in- 
cluding energies in the vicinity of h max where our phe- 
nomenon takes place. Correspondingly, the integral in 
(11) over A acquires the meaning of the integration over 
energy: dX => (dX/dE)dE = (dX/dl) dl. 

Consider (27) in some more details. The transition 
coefficient, \a\, on the real axis is of the form: 



a 



= (1 + e" 



vl/2 



(28) 



It has the branch-cuts at points l m — ±(2m + l)i (m = 
0, 1, . . .) which should be chosen along imaginary axis, 
for \a\ to be the single valued analytical function inside a 
strip near the real axis. The phase, 6(1), displays similar 
properties at |/| ~ 1 (see S(l) (31)). As for p(l) (26), the 
square root in (26) has singularities at: 



l' m = ±(2m + l)i - - ln(sin 2 2ttk) 

7T 



(29) 



Bending the contour of integration over I into the upper 
half- plane, the main contribution into M osc (11) is to be 
due to "nearest" singularities, at l m and l' m . 

These values have the scale |A| ~ 1, which, being re- 
written in terms of the energy, (24) , results in: 



\s\~Lo c a}t 2 ~MW^) 1/2 <ft 



(30) 



Eq. (30) justifies our implicit assumptions above that the 
effect is caused by the main peak, h max . Other peaks, 
although being of the same scale, h ~ VF/d, do not con- 
tribute to the effect produced by the fine structure (30) 
of the levels developing in the very vicinity of h max . 

The subsequent analysis although simple and straight- 
forward, involves somewhat lengthy mathematical details 
to be published elsewhere. Below we provide only a 
sketch of a few major steps, before discussing the result. 



Exact eq. (27), has no non-analytic Hnl-term in 
(19) at small I. At I <~ 1 one gets (c <~ 1): 



S(0 = 5(27r,0) + - 




(31) 



The important step is to calculate |it/| 2 in (11). In 
particular, the normalization is to be known with bet- 
ter accuracy than provided by the WKB-expressions of 
eq. (16, 16'). This can be achieved in frameworks of the 
above method. Fortunately however, properties of the 
Bloch functions for a one-dimensional periodical poten- 
tial are well studied. With the help of eq. (4.18) in ref. 
[8] and our eqs. (21) we derived: 



ui 



| 2 = - — ira 1 / 2 sin 2ttk 



dR 



(32) 



where R(l) is the R.H.S. of eq. (23). Expression for 
M osc (11) is given as a sum of integrals, Ik, which may 
be re- written as integrals over the energy variable, i. 



/+oo 
-OO 



dl 



(hi 2 ) 



dl 



(33) 



The constant, i.e. 1/2, having been eliminated in (33), 
the path of integration over real I can be transformed 
into two contours, C\ and C*2, encircling the two branch- 
cuts caused by singular points of l m , and l' m of 
eq. (28, 29). With integrations over C\ and C2 now 
running parallel to the imaginary axis integrals become 
rapidly (exponentially) convergent, and integration back 
by parts is allowed. One obtains two integrals of the 
form: 



iK-ira 1 / 2 sin 2ttk 



LK 



ex.piK®(l)dl 
(sin 2 2ttk + e -^)V2 



(34) 



With <&(£) given by (27) and (31), and the proper def- 
inition of the branches on the cuts, integrations can 
be completed analytically, particularly so in the limit 

m (\/W?o) 3> 1. Without providing here the result ex- 
plicitly, we use this latter assumption to discuss the order 
of magnitude for Ik=i- 



Ik=\ oc a 



1/2 iS(2irfi) 



(Wd) 



1/4 



AA, c (Wd) 7/ V s ^°> 



(35) 



The numerical coefficient in (35) is large but there is a 
phase factor containing 



(-271 



S(2tt,0) = v^AK / {hmax ~ h{ip))^ 2 dip (36) 
Jo 

and the result has to be averaged yet over all trajectories. 



Our estimates proceed in the following way. Since the 
Larmor orbit is so large (r^ 3> d), the maximum in h{(p) 
must be practically the same for any trajectory, while 
its position, ip , varies. The position, <p , played no role 
in our analysis above. The integral in (36) is positive, 
therefore we write 



iS'(27r,0) 



where 5(2^,0) is an average, i.e. just a phase (al- 
though unknown), while SS fluctuates around 5 f (27r,0). 

Since 5(2tt,0) ~ A/lu c {Co/d) 1/2 is large, (e iSS ) can be 
calc ul ated usin g for SS the Gaussian distribution with 



SS 2 = S(2ir,0). For M^JM^ it gives: 

AK (£ /d) 7/4 exp \-A/u c (^/d) 1/2 



(37) 



The exponent (37) remains, of course, yet too small 
to expect that the dHvA effect may be measured ex- 
perimentally in the regime (1). Nevertheless, we have 
demonstrated the effect of the local density of states char- 
acteristic of the superconducting state: an effective Din- 
gle's temperature becomes reduced as A =^ A^o/d) 1 ^ 2 
even for the isotropic BCS-like superconductor. In the 
anisotropic case such a reduction is expected to be 
stronger [9] A A(£o/d), in the same regime (1). 

To conclude, we demonstrated that, unlike normal 
metals where the dHvA-effect is caused by electron lev- 
els crossing the chemical potential, in the superconduct- 
ing heterogeneous state oscillations are caused by lev- 
els crossing the energy threshold separating "localized" 
and itinerant states. We argue that this mechanism may 
be responsible for pronounced oscillations in the regime 
B ~ H C 2 where not only is the gap weaker, but is het- 
erogeneous (periodic) itself, together with distributions 
of the field and currents. 
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